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The coefficients of different combinations of terms of the liquid drop model have been determined 
by a least square fitting procedure to the experimental atomic masses. The nuclear masses can also 
be reproduced using a Coulomb radius taking into account the increase of the ratio R0/A 1 / 3 with 
increasing mass, the fitted surface energy coefficient remaining around 18 MeV. 
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To predict the stability of new nuclides both in the 
superheavy element region and the regions close to the 
proton and neutron drip lines continuous efforts are still 
needed to determine the nuclear masses and therefore 
the binding energies of such exotic nuclei. Within a 
modelling of the nucleus by a charged liquid drop, semi- 
macroscopic models including a pairing energy have been 
firstly developed to reproduce the experimental nuclear 
masses 0, Q • The coefficients of the Bethe-Weizsacker 
mass formula have been determined once again recently 
0- To reproduce the non-smooth behaviour of the 
masses (due to the magic number proximity, parity of 
the proton and neutron numbers,...) and other micro- 
scopic properties, macroscopic-microscopic approaches 
have been formulated, mainly the finite-range liquid drop 
model and the finite-range droplet model [4,]. Nuclear 
masses have also been obtained accurately within the 
statistical Thomas-Fermi model with a well-chosen ef- 
fective interaction H |(| . Microscopic Hartree-Fock self- 
consistent calculations using mean-fields and Skyrme or 
Gogny forces and pairing correlations 0,0 as well as rel- 
ativistic mean field theories || have also been developed 
to describe these nuclear masses. Finally, nuclear mass 
systematics using neural networks have been undertaken 
recently |To| . 

The nuclear binding energy B„ UC ;(A,Z) which is the 
energy necessary for separating all the nucleons consti- 
tuting a nucleus is connected to the nuclear mass M n m 
by 

B nud (A, Z) = Zm P + Nm N - M n . m (A, Z). (1) 

This quantity may thus be easily derived from the exper- 
imental atomic masses as published in EH since : 

M n . m {A, Z) = M a . m {A, Z) - Zm e + B e {Z) (2) 

while the binding energy B e (Z) of all removed electrons 
is given by [T^ 



B e (Z) = a e iZ^ + b e iZ bdb , 



(3) 



witha e/ = 1.44381 x 1(T 5 MeV and b el = 1.55468 x lO" 12 
MeV. 

The fission, fusion, cluster and a decay potential barri- 
ers are governed by the evolution of the nuclear binding 
energy with deformation. It has been shown that four 
basic terms are sufficient to describe the main features 
of these barriers [3 El El El El El : the volume, 
surface, Coulomb and nuclear proximity energy terms 
while the introduction of the shell and pairing energy 
terms is needed to explain structure effects and improve 
quantitatively the results. Other terms have been pro- 
posed to determine accurately the binding energy and 
other nuclear characteristics : the curvature, A , proton 
form factor correction, Wigner, Coulomb exchange cor- 
rection,... energy terms 0. 

The purpose of the present work is to determine the 
coefficients of different combinations of terms of the liq- 
uid drop model by a least square fitting procedure to the 
experimentally available atomic masses |ll| and to study 
whether nuclear masses can also be reproduced using, for 
the Coulomb energy, a radius which takes into account 
the small decrease of the density with increasing mass 
and to determine the associated surface energy coeffi- 
cient. The theoretical shell effects given by the Thomas- 
Fermi model (7 th column of the table in [5j and have 
been selected since they reproduce nicely the mass decre- 
ments from fermium to Z = 112 |l9j . They are based on 
the Strutinsky shell-correction method and given for the 
most stable nuclei in the appendix. The masses of the 
1522 nuclei verifying the two following conditions have 
been used : N and Z higher than 7 and the one standard 
deviation uncertainty on the mass lower than 100 keV 

m 

The following expansion of the nuclear binding energy 
has been considered 

B nuc i =a v (l- k v I 2 )A - a s (l - k s I 2 )A 2 / 3 - (4) 
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The nuclear proximity energy term does not appear since 
its effect is effective only for necked shapes but not 
around the ground state. The first term is the volume 
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energy and corresponds to the saturated exchange force 
and infinite nuclear matter. In this form it includes the 
asymmetry energy term of the Bethe-Weizsacker mass 
formula via the relative neutron excess I=(N-Z)/A. The 
second term is the surface energy term. It takes into 
account the deficit of binding energy of the nucleons at 
the nuclear surface and corresponds to semi-infinite nu- 
clear matter. The dependence on / is not considered in 
the Bethe-Weizsacker mass formula. The third term is 
the Coulomb energy. It gives the loss of binding energy 
due to the repulsion between the protons. In the Bethe- 
Weizsacker mass formula the proportionality to Z(Z-l) is 
assumed. 

The pairing energy has been calculated using 

Ep a i r — ~a p / A 1 / 2 for odd Z, odd N nuclei, 

Epair — for odd A, (5) 
Epair — cip/A 1 / 2 for even Z, even N nuclei. 

The a p = 11 value has been adopted following first fits. 
Other more sophisticated expressions exist for the pairing 
energy 0,0. 

The sign for the shell energy term comes from the 
adopted definition in 0. It gives a contribution of 
12.84 MeV to the binding energy for 208 Pb for exam- 
ple. The curvature energy afcA 1 / 3 is a correction to the 
surface energy appearing when the surface energy is con- 
sidered as a function of local properties of the surface and 
consequently depends on the mean local curvature. The 
aoA term appears when the surface term of the liquid 
drop model is extended to include higher order terms in 
A^ 1 / 3 and /. The last but one term is a proton form- 
factor correction to the Coulomb energy which takes into 
account the finite size of the protons. The last term is 
the Wigner energy 4, 20] which appears in the counting 
of identical pairs in a nucleus, furthermore it is clearly 
called for by the experimental masses. 

In Table 1 the improvement of the experimental data 
reproduction when additional terms are added to the 
three basic volume, surface and Coulomb energy terms 
is displayed when the nuclear radius is calculated by the 
formula Ro=roA 1 / 3 . 

The root-mean-square deviation, defined by : 

ff2 = S [M Th - M Exp } 2 
n 

has been used to compare the efficiency of these different 
selected sets of terms. 

To follow the non smooth variation of the nuclear 
masses with A and Z the introduction of the shell en- 
ergy is obviously needed as well as that of the pairing 
term though its effect is smaller. The curvature term 
and the constant term taken alone do not allow to better 
fit the experimental masses while the proton form fac- 
tor correction to the Coulomb energy and the Wigner 
term have separately a strong effect in the decrease of 
a. When both these two last contributions to the bind- 
ing energy are taken into account a = 0.58 MeV which 



is a very satisfactory value 0, 0, |2(j ■ The addition to 
the proton form factor and Wigner energy terms of the 
curvature energy or constant terms taken separately or 
together does not allow to improve a. Furthermore, a 
progressive convergence of a& and ao is not obtained and 
strange surface energy coefficients appear. Due to this 
strong variation and lack of stability of the curvature 
and constant coefficient values it seems preferable to ne- 
glect these terms since the accuracy is already correct 
without them. Disregarding the last line a good stability 
of the volume & v and asymmetry volume k v constants is 
observed. The variation of the surface coefficient is larger 
but a s reaches a maximum of only 18 MeV. As it is well 
known the surface asymmetry coefficient k s is less easy 
to precise. Invariably ro has a value of around 1.22 fm 
within this approach. 

For the Bethe-Weizsacker formula the fitting procedure 
leads to 

B nucl (A, Z) = 15.7827A - 17.9042A 2/3 (7) 
-0.72404 ^^3 1} - 23.7193/ 2 A + E pair - E shell 

with cr=1.17 MeV. That leads to r =1.193 fm and 
k t , = 1.505. The non dependence of the surface energy 
term on the relative neutron excess I explains the a value. 

The formula Ro=r A 1 / 3 does not reproduce the small 
decrease of the density with increasing mass [2(J . In pre- 
vious works HUll 111 13 13 111 the formula 

R = l.28A 1/3 -Q.76 + 0.8A- 1/3 (8) 

proposed in Ref. [2l|] for the effective sharp radius has 
been retained to describe the main properties of the fu- 
sion, fission, cluster and a emission potential barriers in 
the quasi-molecular shape path. It leads for example to 
r =1.13 fm for 48 Ca and r =1.18 fm for 248 Cm. The 
fit of the nuclear binding energy by the expression (6) is 
reconsidered in the table 2 in calculating the Coulomb en- 
ergy by the formula f 128A i/3_o.76+o.8A-i/3 witnout ad- 
justable parameter. The behaviour of the different com- 
binations of terms is about the same as in table 1. It 
seems also preferable to disregard the curvature and con- 
stant contributions. Then a takes the value 0.60 MeV 
to be compared to 0.58 MeV in the first table when the 
Coulomb energy contains an additional adjustable pa- 
rameter. 

The figure shows that the difference between the theo- 
retical and experimental masses never exceeds 2.25 MeV 
and is less than 1.15 MeV when A is higher than 100. In 
the Generalized Liquid Drop Model 0, 0] the selected 
values are a„=15.494 MeV, a s =17.9439 MeV, k t ,=1.8 and 
k s =2.6. They are close to the ones given <r=0.60 MeV in 
table 2. 

As a conclusion, the most important result is that 
it is possible to reproduce the nuclear masses in taken 
a realistic formula for the effective sharp radius given 
Ro/A 1 / 3 ^!.! fm for the lightest nuclei and 1.18 fm for 
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TABLE I: Dependence of the energy coefficient values (in MeV or fm) on the selected term set including or not the pairing and 
theoretical shell energies and root mean square deviation. The Coulomb energy is determined by a c -^-^ with a c = 3e 2 /5ro. 
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1.7013 
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1.3175 


1.2046 
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15.4496 


1.8376 


17.3525 


2.0990 


1.2249 






-0.93794 


27.2843 


y 
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15.5518 


1.8452 
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2.0910 


1.2176 


-1.1207 




-0.97833 


28.4963 
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1.8462 
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28.6061 


y 


y 
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1.8373 


11.2996 


2.7852 


1.2505 


20.9967 


-25.7261 


- 0.6862 


26.4072 


y 


y 


0.57 



TABLE II: Dependence of the energy coefficient values (in MeV) on the selected term set including or not the pairing and shell 
energies and the corresponding root mean square deviations. The Coulomb energy coefficient is not adjusted and is determined 



u v 3 e'Z" 

y 5 1.28A 1 /3-o.76+0.8A- 1 /^ " 
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k s 
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15.9622 
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18.0108 


1.0627 
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2.3757 
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31.0077 
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-1.71391 
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y 


0.58 




Wigner energy terms. 

APPENDIX A: TABLE III 
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FIG. 1: Difference (in MeV) between the theoretical and ex- 
perimental masses for the 1522 nuclei as a function of the 
mass number. 



the heaviest ones while the surface energy coefficient re- 
mains around 18 MeV and the surface-asymmetry co- 
efficient around 2.5. Besides the values of the volume, 
surface and Coulomb energy coefficients the accuracy of 
the fitting depends also strongly on the selected shell and 
pairing energies as well as on the proton form factor and 
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TABLE III: Theoretical shell energy (in MeV) extracted from 5] (7 th column) at the ground state of nuclei for which the 
half-life is higher than 1 ky. 



w B 


1X B 


12 C 


13 c 


14 C 


14 7V 


15 N 


16 G 


17 o 


18 


19 p 


20 iV e 


21 Ne 


22 iVe 


2.36 


0.74 


-0.69 


-1.03 


-0.56 


-1.33 


-0.87 


-0.45 


1.19 


1.3 


2.76 


2.81 


2.82 


2.19 


23 TV a 


2A Mg 


25 Mg 


26 Mg 


26 Al 


27 Al 


28 Si 


29 Si 


30 Si 


31p 


32 S 


33 S 


34 5 


36 S 


2.22 


1.64 


1.77 


0.64 


1.89 


0.79 


-0.26 


-0.25 


0.22 


0.22 


0.66 


0.86 


1.13 


1.22 


35 CI 


36 CI 


37 Cl 


36 Ar 


38 Ar 


40 Ar 


39 K 


40^ 


41 K 


40 Ca 


41 Ca 


42 Ca 


43 Ca 


44 Ca 


1.32 


1.48 


1.40 


1.57 


1.64 


2.45 


1.79 


2.58 


2.58 


1.71 


2.49 


2.49 


2.16 


1.7 


46 Ca 


48 Ca 


45 Sc 


46 Ti 


47 Ti 


48 T - 


49 Ti 


50 T . 


50 y 


51 y 


50 Cr 


52 Cr 


53 Cr 


54 Cr 


0.69 


-0.82 


2.43 


2.44 


1.95 


1.44 


0.81 


-0.05 


0.54 


-0.31 


0.74 


-0.69 


-0.38 


0.74 


53 Mn 


55 Mn 


54 Fe 


56 Fe 


57 Fe 


58 Fe 


60 Fe 


59 Co 


58 Ni 


59 Ni 


m Ni 


61 Ni 


a2 Ni 


64 /Vi 


-1.11 


0.68 


-1.54 


0.05 


0.72 


1.22 


2.07 


0.77 


-1.58 


-0.68 


-0.16 


0.59 


1.1 


1.63 


63 Cu 


65 Cu 


64 Zn 


66 Zn 


m Zn 


68 Zn 


70 Zn 


69 Ga 


71 Ga 


70 Ge 


72 Ge 


73 Ge 


74 Ge 


76 Ge 


1.86 


2.33 


2.53 


2.89 


3.16 


2.99 


2.94 


3.79 


3.71 


4.13 


4.08 


4.19 


3.82 


2.53 


75 As 


74 Se 


76 se 


77 Se 


78 Se 


79 Se 


80 Se 


82 Se 


79 Br 


81 Br 


so Kr 


81 Kr 


82 Kr 


83 Kr 


4.08 


4.42 


4.08 


4.06 


3.27 


2.87 


1.89 


0.38 


4.07 


2.28 


4.39 


3.77 


2.74 


1.65 


84 A> 


86 Kr 


85 Kb 


87 Rb 


86 5r 


87 Sr 


88 Sr 


89y 


90 Zr 


91 Zr 


92 Zr 


93 Zr 


94 Zr 


96 Zr 


0.96 


-0.40 


1.13 


-0.35 


0.79 


0.05 


-0.97 


-1.19 


-1.63 


-0.47 


0.46 


1.53 


2.54 


3.49 


92 Nb 


93 Nb 


94 Nb 


92 Mo 


93 Mo 


94 Mo 


95 Mo 


96 Mo 


97 Mo 


98 Mo 


100 Mo 


97 Tc 


98 Tc 


"Tc 


-0.57 


0.44 


1.51 


-2.12 


-1.07 


-0.12 


0.97 


1.79 


2.46 


2.98 


3.62 


1.26 


2.05 


2.64 


96 Ru 


98 Ru 


"Ru 


100 Ru 


101 Ru 


102 D 

Ru 


104 r> 

Ru 


103 Rh 


102 Pd 


104 Pd 


105 Pd 


106 Pd 


107 Pd 


108 Pd 


-1.11 


0.57 


1.36 


2.00 


2.54 


2.98 


3.49 


2.44 


0.63 


1.83 


2.39 


2.80 


3.08 


3.34 


110 Pd 


107 Ag 


109 Ag 


106 Cd 


108 Cd 


110 Cd 




112 Cd 


113 Cd 


114 Cd 


116 Cd 


113 /n 


115 /n 


112 Sn 


3.42 


2.20 


2.91 


0.31 


1.35 


2.11 


2.42 


2.52 


2.61 


2.50 


2.26 


1.83 


1.97 


0.37 


114 Sn 


115 Sn 


116 Sn 


117 Sn 


118 Sn 


119 Sn 


120 Sn 


122 Sn 


124 Sn 


12a Sn 


121 56 


123 S6 


120 T( , 


122 Te 


0.81 


1.03 


0.94 


0.96 


0.75 


0.70 


0.19 


-0.99 


-2.51 


-4.36 


0.76 


-0.16 


2.08 


1.55 


123 Te 


124 Te 


125 Te 


126 Te 


128 Te 


130 Te 


127 j- 


129 j 


128 Xe 


129 Xe 


130 Xe 


131 Xe 


132 Xe 


134 Xe 


1.30 


0.66 


0.25 


-0.62 


-2.46 


-4.74 


0.35 


-1.29 


1.01 


0.48 


-0.31 


-1.11 


-2.36 


-4.92 


136 Xe 


133 Cs 


135 Cs 


132 Ba 


134 Ba 


135 Ba 


136 Ba 


137 Ba 


138 Ba 


137 La 


138 La 


139 La 


136 Ce 


138 Ce 


-7.2 


-1.28 


-3.90 


0.93 


-0.55 


-1.45 


-3.01 


-4.18 


-5.29 


-2.22 


-3.37 


-4.50 


0.70 


-1.57 


140 Ce 


142 Ce 


141 Pt 


142 Nd 


143 Nd 


144 Nd 


145 Nd 


146 Nd 


148 7Vd 


150 7Vd 


144 Sm 


146 Sm 


147 Sm 


148 Sm 


-3.86 


-2.06 


-3.26 


-2.88 


-2.14 


-1.04 


0.10 


0.56 


0.85 


0.54 


-2.28 


-0.41 


0.67 


1.12 


149 Sm 


150 5m 


152 Sm 


154 5m 


151 Eu 


153 Eu 


150 Gd 


152 Gd 


154 Gd 


155 Gd 


156 Gd 


157 Gd 


158 Gd 


160 Gd 


1.22 


1.31 


0.90 


0.38 


1.39 


1.02 


1.30 


1.59 


1.33 


1.05 


0.89 


0.62 


0.56 


0.21 


159 Tb 


154 Dy 


156 Dy 


158 Dy 


160 Dy 


1G1 Dy 


162 Dy 


163 Dy 


W Dy 


165 Dy 


163 Ho 


1(i5 Ho 


162 Er 


164 Er 


0.64 


1.63 


1.56 


1.24 


0.92 


0.64 


0.47 


0.16 


-0.06 


-0.41 


0.46 


-0.12 


1.20 


0.70 


lm Er 


167 Er 


168 Er 


170 Er 


169 Tm 


168 Yb 


170y b 


171 Yb 


172y b 


173 Yb 


174y 6 


176y b 


175 Lu 


176 Lu 


0.07 


-0.37 


-0.54 


-1.08 


-0.60 


0.32 


-0.34 


-0.76 


-0.94 


-1.32 


-1.30 


-1.74 


-1.23 


-1.62 


174 Hf 


176 Hf 


177 Hf 


178 Hf 


179 Hf 


180 Hf 


182 Hf 


181 Ta 


180 W 


182 w 




184 W 


186 W 


185 Re 


-0.38 


-0.90 


-1.33 


-1.53 


-1.97 


-1.99 


-2.16 


-2.02 


-1.21 


-1.71 


-2.00 


-2.02 


-2.38 


-2.19 


187 7-» 

Re 


184 /■ \ 

Os 


186 ^ 

Os 


187 ^ 

Os 


188 s-\ 


189 

Os 


190 ^ i 

Os 


192 

Os 


191 T 

• Ir 


193 t 


190 t-> i 


192 7-1/ 


194 Ti i 
M Pt 


195 7-1/ 
J0 Pt 


-2.48 


-1.61 


-1.88 


-2.16 


-2.08 


-2.43 


-2.47 


-3.50 


-2.54 


-3.62 


-0.97 


-2.01 


-3.31 


-4.04 


196 Pt 


198 Pi 


197 Au 


196 Hg 


198 Hg 


™Hg 


200 Hg 


201 Hg 


202 Hg 


204 Hg 


203y; 


205 rp^ 


202 Pb 


204 Pb 


-4.80 


-6.11 


-5.56 


-4.51 


-5.99 


-6.75 


-7.52 


-8.37 


-9.11 


-10.69 


-9.97 


-11.58 


-8.22 


-10.02 


2 05pb 


W6 Pb 


207 Pb 


208p b 


208 Bi 


209 Bi 


226 Ra 


229 Th 


230 Th 


232 Th 


231 Pa 


233 jj 


234 U 


235 jj 


-11.00 


-11.82 


-12.68 


-12.84 


-11.70 


-11.95 


-0.30 


-0.52 


-0.43 


-0.60 


-0.79 


-1.27 


-1.23 


-1.46 


236 jj 


238 


236 Np 


237 Np 


239 Pu 


240 Pu 


242 Pu 


244 Pu 


243 Am 


245 Cm 


246 Cm 


247 Cm 


248 Cm 


247 Bk 


-1.30 


-1.27 


-1.85 


-1.74 


-2.12 


-1.95 


-1.99 


-2.08 


-2.44 


-3.05 


-2.96 


-3.17 


-3.00 


-3.46 
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